Identification of the Molecule-Metal Bonding Geometries of Molecular Nanowires 
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Molecular nanowires in which a single molecule bonds chemically to two metal electrodes and forms 
a stable electrically conducting bridge between them have been studied intensively for more than 
a decade. However the experimental determination of the bonding geometry between the molecule 
and electrodes has remained elusive. Here we demonstrate by means of ab initio calculations that 
inelastic tunneling spectroscopy (IETS) can determine these geometries. We identify the bonding 
geometries at the gold-sulfur interfaces of propanedithiolate molecules bridging gold electrodes that 
give rise to the specific IETS signatures that were observed in recent experiments. 
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Molecular nanowires in which a single organic molecule 
bonds chemically to two metal electrodes and forms 
a stable electrically conducting bridge between them 
have attracted a great deal of attention [1] because of 
their fundamental interest and potential applications as 
single-molecule nanoelectronic devices. However, a single 
molecule located between two electrodes is not accessible 
to scanning microprobes that can measure atomic scale 
structure. Thus definitive determination of the bonding 
geometries at the molecule-electrode interfaces of single- 
molecule nanowires continues to be elusive despite being 
critically important for understanding electrical conduc- 
tion through molecular wires quantitatively, and for gain- 
ing control over their structures for device applications. 

In the case of single-molecule nanowires bridging gold 
electrodes and thiol bonded to them, possible bonding 
geometries include those in which a sulfur atom of the 
molecule is located at a top site over a particular gold 
surface atom or over a bridge site between two gold atoms 
or over a hollow site between three gold atoms. [1] How- 
ever, there have been no direct experiments determining 
which (if any) of these possibilities are actually realized. 
For molecules amine-linked to gold electrodes it has been 
argued that top-site bonding is the most probable [2], but 
direct experimental evidence of this has also been lacking. 

The excitation of molecular vibrations (phonons) by 
electrons passing through single-molecule nanowires gives 
rise to conductance steps in the low temperature current- 
voltage characteristics of the nanowires. Inelastic tunnel- 
ing spectroscopy (IETS) experiments have detected these 
steps and measured the energies of the emitted phonons 
[3]. These experiments proved that particular molecu- 
lar species are involved in electrical conduction through 
metal-molecule-metal junctions. Density functional the- 
ory (DFT) based simulations [4-6] have accounted for 
many features of the IETS data. However, the possibil- 
ity that IETS might identify the bonding geometries at 
the molecule-metal interfaces and thus resolve the long 
standing problem of the molecule-electrode structure has 
not been investigated in the literature. We explore it in 
this paper and identify for the first time the molecule- 
contact bonding geometries that were realized in exper- 
iments on an organic molecule bridging metal contacts. 



We consider one of the simplest organic molecules, 1,3- 
propanedithiolate (PDT), bridging gold electrodes for 
which detailed experimental IETS data is available. [7] 

We focus on inelastic tunneling processes that are sen- 
sitive to the structure of the gold-sulfur interfaces, i.e., 
those that involve excitation of vibrational modes with 
strong amplitudes on the sulfur atoms. Therefore it is 
necessary to calculate accurate equilibrium geometries 
and also the frequencies and atomic displacements from 
equilibrium for the vibrational modes of the whole sys- 
tem, including both the molecule and the gold electrodes. 
We do this by performing ab initio DFT calculations [8] 
for extended molecules consisting of a PDT molecule and 
two finite clusters of gold atoms that the molecule con- 
nects, relaxing this entire structure. By carrying out 
systematic DFT studies of extended molecules with gold 
clusters of different sizes (up to 13 gold atoms per cluster) 
we establish that our conclusions arc independent of the 
cluster size for the larger clusters that we study and thus 
are applicable to molecules bridging the nanoscalc tips of 
experimentally realized macroscopic gold electrodes. 

We found extended molecules whose sulfur atoms bond 
to the gold clusters over bridge sites to have lower ener- 
gies than those with top site bonding by at least 0.76 
cV. Extended molecules for which DFT geometry relax- 
ations were started with the sulfur atoms over hollow 
sites on the surfaces of close packed gold clusters invari- 
ably relaxed to bridge bonding site geometries. While 
we found it possible to generate examples of relaxed ex- 
tended molecules with each sulfur atom bonding to three 
gold atoms (i.e., hollow site bonding) the structures of 
the gold clusters near these bonding sites were much more 
open than that of a fee gold crystal. The energies of these 
structures were also higher than those of either bridge or 
top site bonded extended molecules. Because of the much 
greater fragility and higher energies of hollow site bonded 
structures relative to bridge and top site bonding we will 
focus primarily on bridge and top site bonding, but will 
revisit hollow site bonding near the end of this paper. 

After identifying the normal modes of the extended 
molecule that have the largest vibrational amplitudes on 
the sulfur atoms and calculating the frequencies of those 
modes as is outlined above, we determine which of these 
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modes has the strongest IETS intensity, i.e, gives rise to 
the largest conductance step height as the bias voltage 
applied across the extended molecule is varied. 

We calculate the IETS intensities perturbatively in the 
spirit of an approach proposed by Troisi et al.[9] who 
transformed the problem of calculating IETS intensities 
into an elastic scattering problem. However, we formu- 
late our IETS intensities explicitly in terms of elastic elec- 
tron transmission amplitudes tj\ through the molecular 
wire. We find the height 5g a of the conductance step due 
to the emission of phonons of vibrational mode a to be 
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at low temperatures. Here ^'({0}) is the elastic elec- 
tron transmission amplitude through the molecular wire 
in its equilibrium geometry from state i with velocity 
Vi in the electron source to state j with velocity Vj in 
the electron drain. d„ Q represents the displacements 
from their equilibrium positions of the atoms n of the 
extended molecule in normal mode a normalized so that 
Yl n m ndna' '^na = <5a',a where m„ is the mass of atom n. 
uj a is the frequency of mode a. tjl({Ad na }) is the elastic 
electron transmission amplitude through the molecular 
wire with each atom n displaced from its equilibrium po- 
sition by Ad na where A is a small parameter. Our formal 
derivation of Eq. (1) will be presented elsewhere; here we 
point out that Eq. (1) is plausible intuitively in a sim- 
ilar way to the result of Troisi et al.[9] : Eq. (1) states 
that in the leading order of perturbation theory the scat- 
tering amplitude for inelastic transmission of an electron 
through the molecular wire is proportional to the change 
in the elastic amplitude for transmission through the wire 
if its atoms are displaced from their equilibrium positions 
as they are when vibrational mode a is excited. 

In our evaluation of tj\ in Eq. (1) the coupling of the 
extended molecule to the macroscopic electron reservoirs 
was treated as in previous work[10-12] by attaching a 
large number of semi-infinite quasi-one-dimensional ideal 
leads to the valence orbitals of the outer gold atoms of the 
extended molecule. The transmission amplitudes i|' were 
then found by solving the Lippmann-Schwinger equation 
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where |$o) is an electron eigenstate of an ideal left lead 
that is decoupled from the extended molecule, Go{E) 
is the Green's function for the decoupled system of the 
ideal leads and the extended molecule, W is the cou- 
pling between the extended molecule and leads, and 
'I') is the scattering eigenstate of the complete coupled 
system. In evaluating Gq(E) semi-empirical extended 
Hiickcl theory[13, 14] was used to model the electronic 
structure of the extended molecule. Previous calcula- 
tions based on extended Hiickel theory have yielded elas- 
tic tunneling conductances in agreement with experiment 
for molecules thiol bonded to gold electrodes [12, 15, 16] 
and have accounted for transport phenomena observed 



in molecular arrays on silicon [10] as well as electrolumi- 
nescence data[17], current-voltage characteristics[17] and 
STM images [18] of molecules on complex substrates. 

We calculated the zero bias tunneling conductances for 
gold-PDT-gold molecular wires from the Landauer for- 
mula g = goYlij \ t fi({°})\ 2v ]/ v i with 5o = 2e 2 //i, evalu- 
ating the elastic transmission amplitudes i|'({0}) as de- 
scribed above and found values in the range g =0.0012- 
0.0014 5o and 0.0012-0.00205 for top site and bridge site 
bonded molecules, respectively. The degree of agree- 
ment between these theoretical values and the experi- 
mental value[7] 0.006 ± 0.002t/o is typical of that in the 
literature[19] on molecules thiol bonded to gold elec- 
trodes; as in previous studies [19] comparing the clastic 
conductance calculations with experiment does not reveal 
which bonding geometry was realized experimentally. 

We then calculated the vibrational modes, their fre- 
quencies and IETS intensities as described above for 
many gold-PDT-gold extended molecules with various 
gold-sulfur bonding geometries and compared our results 
with the experimental data of Hihath et al. [7] . Our cal- 
culations showed that the modes with strong amplitudes 
of vibration on the sulfur atoms that have the largest 
IETS intensities fall within the phonon energy range of 
a prominent feature of the experimental IETS phonon 
histogram[7] that extends from 39 to 52 meV. 

We show the vibrational modes that we find in this 
energy range for representative top site and bridge site 
bonding geometries in Fig. 1. The calculated IETS spec- 
tra (phonon energies and IETS intensities) in the same 
energy range are shown in Fig. 2 for several extended 
molecules with gold clusters of various sizes together with 
the experimental IETS phonon mode histogram. [7] 

The mode with the strongest calculated IETS intensi- 
ties in Fig. 2 is mode I shown in the top row of Fig. 1. In 
this mode the sulfur atoms have the strongest vibrational 
amplitudes and move in antiphase, approximately along 
the axis of the molecule. The other mode in this energy 
range is mode II shown in the lower row of Fig. 1. It is 
similar to mode I except that the sulfur atoms move in 
phase with each other. As is seen in Fig. 2 the calculated 
IETS intensities for mode II are much weaker than those 
for mode I for both top and bridge site bonding. The 
reason for this difference between the IETS intensities 
of modes I and II can be understood intuitively by con- 
sidering the nature of the motion in relation to Eq. (1): 
Since in mode I the two sulfur atoms move in antiphase 
the gold-sulfur distances for both sulfur atoms either in- 
crease or decrease together as the extended molecule vi- 
brates. These distances can be regarded as the widths of 
tunnel barriers between the molecule and the two gold 
electrodes. Thus the motions of the two sulfur atoms 
act in concert to widen or narrow both tunnel barriers 
together and therefore to weaken or strengthen the elec- 
tron transmission amplitude through the molecular wire. 
Thus the magnitude of the difference between the elas- 
tic electron transmission amplitudes through the molec- 
ular wire in its equilibrium and vibrating geometries 
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FIG. 1. Calculated vibrational modes in phonon energy range 
from 39 to 52 meV for trans-PDT bridging gold nano-clusters 
with sulfur atoms bonded to gold in top-site and bridge-site 
geometries. Red arrows show (arbitrarily normalized) atomic 
displacements. Mode I has the stronger IETS intensity. [20] 
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FIG. 2. Calculated phonon energies and IETS intensities for 
trans-PDT molecules linking pairs of gold clusters with be- 
tween 10 and 13 Au atoms in each cluster. Results are shown 
for both sulfur atoms bonding to the gold in top-site and 
bridge-site geometries and for top site bonding to one gold 
cluster and bridge site bonding to the other. The two ellipses 
enclose the calculated type I mode IETS spectra for pure 
bridge and mixed top-bridge bonding geometries respectively 
for extended molecules with gold clusters containing various 
numbers of gold atoms. The experimental IETS phonon mode 
histogram of Hihath et al. [7] is shown in (darker, lighter) grey 
for (positive, negative) bias voltages. 



\tfidAdna}) - £|!({0})| in Eq. (1) is enhanced. By con- 
trast in mode II when the gold-sulfur distance for one sul- 
fur atom increases the gold-sulfur distance for the other 
sulfur atom decreases. Thus the effects of the motions 
of the two sulfur atoms on the elastic transmission am- 
plitude through the molecular wire tend to cancel. Thus 
Itf^Adna}) - ^-({0})| in Eq. (1) is smaller for mode 



for 



II than mode I and therefore the IETS intensity 6g, 
mode II is much smaller as is seen in Fig. 2. 

Comparison of the calculated IETS spectra in Fig. 2 
with the experimental phonon mode histogram [7] indi- 
cates that vibrational mode I contributed most of the 
counts recorded in the histogram in the energy range 
shown. This is consistent with the fact that the calcu- 
lated IETS intensities for mode I are much stronger than 



those for mode II and therefore mode I should be more 
readily detected in experimental IETS measurements. 

The theoretical results for the dominant IETS mode 
I in Fig. 2 reveal that the prominent features of the 
histogram [7] can be explained as arising from PDT 
molecules that bond to the gold electrodes in different 
ways: The main peak in the experimental histogram[7] 
that is centered at ~46 meV matches our theoretical re- 
sult for trans-PDT molecules that bond to both gold elec- 
trodes in the top site geometry. The weaker peak cen- 
tered near 42 meV matches our predictions for molecules 
that bond to one gold electrode in the top site geometry 
and to the other electrode in the bridge site geometry. 
The shoulder of the experimental histogram [7] at lower 
phonon energies between 39.5 meV and 41.5 meV cor- 
responds to our results for molecules bonding to both 
electrodes in the bridge site geometry. 

Notice that our theoretical results for phonon mode I of 
PDT molecules bonded to both electrodes in the top site 
geometry (the feature near 45.5 meV in the theoretical 
spectra in Fig. 2) are very well converged with respect to 
increasing gold cluster size: Both the calculated phonon 
energies and the IETS intensities are almost independent 
of the gold cluster size in the size range shown (10-13 
gold atoms per cluster). The calculated phonon energy 
of this dominant mode matches the phonon energy of the 
main peak of the experimental phonon mode histogram 
in Fig. 2 very well. It is also well separated from the 
calculated phonon energies of the dominant mode I for 
the pure bridge and mixed bridge-top site bonded ge- 
ometries. As a check we calculated the vibrational mode 
energies and IETS intensities for a few examples of ex- 
tended molecules using a different density functional [8] 
and found similar results. Thus our results identify un- 
ambiguously those specific realizations of the molecular 
wire that gave rise to the counts within the main peak 
of the experimental histogram[7] reproduced in Fig. 2 
as being those in which both sulfur atoms bonded to 
the gold electrodes in the top site geometry. The calcu- 
lated mode I phonon energies in Fig. 2 for pure bridge 
site-bonded and mixed bridge-top site-bonded molecular 
wires show more variation with gold cluster size than do 
the calculated mode I phonon energies for molecules in 
the pure top site-bonded geometry. However, the ranges 
in which the calculated energies of the mode I phonons 
for pure bridge site-bonded and mixed bridge/top site- 
bonded molecular wires occur do not overlap. Also, the 
calculated phonon energies for these modes and bonding 
geometries do not show any systematic trend towards 
higher or lower values with increasing cluster size. Thus 
it is plausible that the weaker peak centered near 42 meV 
in the experimental histogram is due to mixed bridge/top 
site-bonded molecular wires and the lower energy shoul- 
der between 39.5 meV and 41.5 meV in the histogram is 
due to pure bridge site-bonded wires. 

As the gold-PDT-gold junction was stretched in the 
experiment [7] the energy of the prominent phonon mode 
in the IETS spectrum was observed to switch from ^42 
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meV to ^46 mcV. It was conjectured [7] that this switch 
may be due to a change in the contact configuration 
or between gauche and trans molecular geometries, but 
no evidence supporting either possibility was offered and 
the contact configurations involved were not identified. [7] 
Here we point out that trans-PDT molecules switch- 
ing from the mixed top-bridge site bonding (calculated 
phonon energy ^42 mcV) to the pure top site bonding ge- 
ometry (calculated phonon energy ~45.5 meV) accounts 
for this observed transition. Our calculated distance be- 
tween the gold clusters of the extended molecule is larger 
for pure top site bonding than for mixed top-bridge bond- 
ing. This is consistent with the transition from mixed 
top-bridge to pure top site bonding occurring as the 
molecular junction is stretched. For longer chain alka- 
nedithiolates, it has been conjectured[21] that switching 
from mixed hollow site-top site bonding, to top site bond- 
ing at both gold electrodes may occur. This conjecture 
does not account for the observed switch from the ^42 
meV mode to the ~46 meV mode in the gold-PDT-gold 
junctions: As noted above, we find hollow-site bonding 
to be much more fragile (and thus less likely to be real- 
ized) than bridge site bonding. We also find the energy 
of the mode of the mixed hollow-top structure with the 
strongest IETS intensity to match neither the ~42 meV 
nor the ^46 meV phonon mode. The theoretical results 
presented above are for trans-PDT molecules. We have 



also studied many gold-PDT-gold structures with molec- 
ular gauche defects and find that such structures also do 
not account for the observed [7] switching from the ~42 
meV phonon mode to the ^46 meV phonon mode. 

In conclusion: We have shown inelastic tunneling spec- 
troscopy to be able to distinguish between different bond- 
ing geometries of the molecule and metal contacts in 
single-molecule molecular wires, an important and previ- 
ously elusive goal in the field of single-molecule nanoelec- 
tronics. We have definitively identified particular realiza- 
tions of gold-propancdithiolatc-gold molecular wires in a 
recent experiment [7] in which the molecule bonded to a 
single gold atom of each electrode. The success of our ap- 
proach rests on the fact that ab initio density functional 
theory calculations of vibrational modes and their fre- 
quencies are known to be accurate because in the Born- 
Oppenheimer approximation they are electronic ground 
state total energy calculations. [1, 5] We rely on trans- 
port calculations only for the identification of the phonon 
mode in a particular frequency range that has the largest 
IETS intensity, and our identification of this mode is also 
supported by physical reasoning. 
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